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1 Introduction 

Soft condensed matter systems are typically multicomponent mixtures of 
macromolecules and simpler components that form complex structures span- 
ning wide ranges of length and time scales [1, 2, 3, 4, 5, 6]. Common classes 
of soft materials are colloidal dispersions [7, 8, 9, 10], polymer solutions 
and melts [11, 12, 13, 14], amphiphilic systems [15, 16], and liquid crys- 
tals [17, 18, 19]. Among these classes are many biologically important systems, 
such as DNA, proteins, and cell membranes. Most soft materials are intrin- 
sically nanostructured, in that at least some components display significant 
variation in structure on length scales of 1-10 nanometers. Many characteristic 
traits of soft condensed matter, e.g., mechanical fragility, sensitivity to exter- 
nal perturbation, and tunable thermal and optical properties, result naturally 
from the mingling of microscopic and mesoscopic constituents. 

The complexity of composition that underlies the rich physical properties 
of soft materials poses a formidable challenge to theoretical and computa- 
tional modelling efforts. Large size and charge asymmetries between macro- 
molecules (e.g., colloidal or polyelectrolyte macroions) and microscopic com- 
ponents (e.g., counterions, monomers, solvent molecules) often render imprac- 
tical the explicit modelling of all degrees of freedom over physically significant 
length and time scales. Model complexity can be greatly reduced, however, by 
pre-averaging (coarse-graining) the degrees of freedom of some of the micro- 
scopic components, thus mapping the original model onto an effective model, 
with a reduced number of components, governed by effective interparticle in- 
teractions. 

The concept of effective interactions has a long history in the statistical me- 
chanics of liquids [20] and other condensed matter systems, dating back over 
a half-century to the McMillan-Mayer theory of solutions [21], the Derjaguin- 
Landau-Verwey-Overbeek (DLVO) theory of charge-stabilised colloids [22, 23], 
and the pseudopotential theory of simple metals and alloys [24, 25]. In mod- 
elling materials properties of simple atomic or molecular liquids and crystals. 
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it is often justifiable to average over electronic degrees of freedom of the con- 
stituent molecules. The coarsci-graincd model then comprises a collection of 
structureless particles interacting via effective intermolecular potentials whose 
parameters depend implicitly on the finer electronic structure. In recent years, 
analogous methods have been carried over and adapted to the realm of macro- 
molecular (soft) materials. 

Effective interparticle interactions prove especially valuable in modelling 
materials properties of soft matter systems, which depend on the collective 
behaviour of many interacting particles. In studies of thermodynamic phase 
behaviour, for example, effective interactions provide essential input to molec- 
ular (e.g., Monte Carlo and molecular dynamics) simulations and statistical 
mechanical theories. While such methods can be directly applied, in princi- 
ple, to an explicit model of the system, brute force applications are, in prac- 
tice, often simply beyond computational reach. Consider, for example, that a 
molecular simulation of only 1000 macroions, each accompanied by as few as 
100 counterions, entails following the motions of 10^ particles, computing at 
each step electrostatic and excludcd-volume interactions among the particles. 
A far more practical strategy applies statistical mechanical methods to an 
effective model of fewer components. 

This chapter reviews the statistical mechanical foundations underlying the- 
ories of effective interparticle interactions in soft matter systems. After first 
identifying and defining the main systems of interest in Sec. 2, several of the 
more common theoretical methods are sketched in Sec. 3. Concise derivations 
are given, in particular, for response theory, density-functional theory, and 
distribution function theory. Although these methods are all well established, 
the interconnections among them are not widely recognized. An effort is made, 
therefore, to demonstrate the underlying unity of these seemingly disparate 
approaches. Practical implementations are illustrated in Sec. 4, where recent 
applications to charged colloids, colloid-polymer mixtures, and polymer solu- 
tions are outlined. In the limited space available, little more than a sample of 
many methods and applications can be included. Complementary perspectives 
and details can be foimd in several excellent reviews [26, 27, 28, 29]. Finally, 
in Sec. 5, a gaze into the (liquid) crystal ball portends an exciting outlook for 
the field. 



2 Systems of Interest 

The main focus of this chapter is effective interactions among macromolecules 
dispersed in simple molecular solvents. The term "simple" here implies sim- 
plicity of molecular structure, not necessarily properties, thereby including 
water - the most ubiquitous, biologically relevant, and anomalous solvent. 
The macromolecules may be colloidal (or nano-) particles, polymers, or am- 
phiphiles and may be electrically neutral or charged, as in the cases of charge- 
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stabilised colloids, polyelectrolytes (including biopolymers), and ionic surfac- 
tants. 

Colloidal suspensions [7] consist of ultra-divided matter dispersed in a 
molecular solvent, and are often classified as lyophobic ("solvent hating") or 
lyophilic ("solvent loving"), according to the case with which the particles 
can be redispersed if dried out. Depending on density, colloidal particles are 
typically nanometers to microns in size - sufficiently large to exhibit random 
Brownian motion, perpetuated by collisions with solvent molecules, yet small 
enough to remain indefinitely suspended against sedimentation. The upper 
size limit can be estimated by comparing the change in gravitational energy 
as a colloid traverses one particle diameter to the typical thermal energy 
ksT at absolute temperature T, where ks is Boltzmann's constant. To better 
appreciate these length scales, consider repeatedly dividing a cube of side 
length 1 cm until reaching first the width of a human hair (10-100 /im), then 
the diameter of a colloid, and finally the diameter of an atom. How many cuts 
are required? 

Polymers [11, 12] are giant chainlikc, branched, or networked molecules, 
consisting of covalently linked repeat units (monomers), which may be all 
alike (homopolymers) or of differing types (heteropolymers) . Polyelectrolytes 
[13, 14] arc polymers that carry ionizable groups. Depending on the nature of 
intramolecular monomer- monomer interactions, polymer and polyelectrolyte 
chains may be stiff or flexible. Flexibility can be quantified by defining a 
persistence length as the correlation length for bond orientations, i.e., the 
distance along the chain over which bond orientations become decorrelated. 
At one extreme, rigid rodlike polymers have a persistence length equal to the 
contour length of the chain. At the opposite extreme, freely-jointed polymers 
are random-walk coils with spatial extent best characterized by the radius of 
gyration, defined as the root-mean-square displacement of monomers from the 
chain's centre of mass. 

Amphiphilic molecules [15, 16] consist of a hydrophilic head group joined to 
a hydrophobic tail group, usually a hydrocarbon chain. The head group may 
be charged (ionic) or neutral (nonionic). When sufficiently concentrated in 
aqueous solution, surfactants and other amphiphiles organize (self-assemble) 
into regular structures that optimize exposure of head groups to the exte- 
rior water phase, while sequestering the hydrophobic tails within. The various 
structures include spherical and cylindrical micelles, bilayers, vesicles (bilayer 
capsules), and microemulsions. Common soap films, for example, are bilayers 
of surfactants (surfactant- water-surfactant sandwiches) immersed in air, while 
biological membranes are bilayers of phospholipids immersed in water. Rela- 
tive stabilities of competing structures are governed largely by concentration 
and geometric packing constraints, as determined by the relative sizes of head 
and tail groups [30]. 

Colloids, polyelectrolytes, and amphiphiles can acquire charge in solution 
through dissociation of ions from chemical groups on the colloidal surfaces, 
polymer backbones, or amphiphile head groups. For ions of sufficiently low va- 
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lence, the entropic gain upon dissociation exceeds the energetic cost of charge 
separation, resulting in a dispersion of charged macroions and an entourage of 
oppositely charged counterions. In an electrolyte solvent, charged macroions 
interact with one another, and with charged surfaces, via electrostatic inter- 
actions that are screened by intervening counterions and salt ions in solution. 
Equilibrium and nonequilibrium distributions of ions are determined by a com- 
petition between entropy and various microscopic interactions [30], including 
repulsive Coulomb and steric interactions and attractive van der Waals (e.g., 
dipole-induced-dipole) interactions. By changing system parameters, such as 
macroion properties (size, charge, composition), and solvent properties (salt 
concentration, pH, dielectric constant, temperature), the range and strength 
of interparticle interactions can be widely tuned. Rational control over the 
enormously rich equilibrium and dynamical properties of macromolecular ma- 
terials relies on a fundamental understanding of the nature and interplay of 
microscopic interactions. 

In all of the systems of interest, the macromolecules possess some degree of 
internal structure. In charged colloids and polyelectrolytes, a multitude of mi- 
croscopic degrees of freedom are associated with the distribution of charge over 
the macroion and the distribution of counterions throughout the solvent. In 
polymer and amphiphilic solutions, the polymer chains or amphiphilic assem- 
blies have conformational freedom. Furthermore, the solvent itself contributes 
a vast number of molecular degrees of freedom. The daunting prospect of ex- 
plicitly modelling multicomponcnt mixtures on a level so fine as to include all 
molecular degrees of freedom motivates the introduction of effective models 
governed by effective interactions. The loss of structural information upon 
coarse graining necessitates an inevitable compromise in accuracy. The art 
of deriving and applying effective interactions lies in crafting approximations 
that are computationally manageable yet capture the essential physics. 

3 Effective Interaction Methods 
3.1 Statistical Mechanical Foundation 

Effective interaction methods have a rigorous foundation in the statistical me- 
chanics of mixtures [31]. These methods rest on the premise that, by averaging 
over the degrees of freedom of some of the components, a multicomponcnt mix- 
ture can be mapped onto an effective model, with a reduced number of compo- 
nents. While the true mixture is subject to bare interparticle interactions, the 
reduced model is governed by coarse-grained, effective interactions. In charged 
colloids, for example, averaging over coordinates of the solvent molecules and 
of the microions (counterions and salt ions) maps the suspension onto an 
effective one-component model of mesoscopic "pseudomacroions" subject to 
microion-induced effective interactions. The bare electrostatic (Coulomb) in- 
teractions between macroions in the true suspension are replaccid by scrcicned- 
Coulomb interactions in the effective model. Similarly, in mixtures of colloids 
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and non- adsorbing polymers, averaging over polymer degrees of freedom leads 
to polymer-induced effective interactions between the colloids. 

Generic Two-Component Model 

We consider here a simple, pairwise-interacting. two-component mixture of Na 
particles of type a and of type b, obeying classical statistics, confined to vol- 
ume V at temperature T. A given system is modelled by three bare interpar- 
ticle pair potentials, Vaa{r), Vbb{r), and Vab{r), assumed here to be isotropic, 
where r is the distance between the particle centres. The model system could 
represent, e.g., a charge-stabilised colloidal suspension, a polyelectrolyte so- 
lution, or various mixtures of colloids, nonadsorbing polymers, nanoparticles, 
or amphiphilic assemblies (micelles, vesicles, etc.). For simplicity, the discus- 
sion is confined to binary mixtures, although the methods discussed below 
easily generalize to multicomponent mixtures. Throughout the derivations, it 
may help to visualize, for concreteness, the a particles as colloids and the b 
particles as counterions. 




Fig. 1. Left: Generic model of a binary mixture of species labelled a and 6. Right: 
Effective one-component model, after coarse-graining of b species, and geometry for 
physical interpretation of response theory. Vectors r and r' define centre-to-centre 
displacements of a particles. Vectors ri, r2, and rs define points at which either the 
"external" potential of the a particles acts or a change is induced in the density of 
b particles (see Sec. 3.2). 

Even before analysing the system in detail, it should be conceptually ap- 
parent that, since particles inevitably influence their environment, the pres- 
ence of particles of one species can affect the manner in which all other parti- 
cles interact. Familiar analogies may be identified in any mixture of interacting 
entities - from inanimate particles to living cells, organisms, and ecosystems. 
In a simple binary mixture, for example, particles of type b can be regarded as 
inducing interactions between a particles. The induced interactions, which act 
in addition to bare aa interactions, depend on both the ab interactions and the 
distribution of b particles. The effective interactions, which are simply sums 
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of induced and bare interactions, may be many-body in character, even if all 
bare interactions arc strictly pairwise, and may depend on the thermodynamic 
state (density, temperature, etc.) of the system. 

Reduction to Effective One-Component Model 

These qualitative observations are now quantified by developing a statistical 
description of the system. We start from the Hamiltonian function H, which 
governs all equilibrium and dynamical properties of the system, and assume 
pairwise bare interactions. The Hamiltonian naturally separates, according to 
H = K + Haa + Hhb + Hab, into the kinetic energy K and three interaction 
terms: 

-ffaa = 2 X! '"aainj), a = a,b (1) 

and 

Hab = J212^''''^^'^^' (2) 
i=l j=l 

where = jr^ — rj| denotes the separation between the centres of particles 
i and j, at positions rj and r^. Within the canonical ensemble (constant Na, 
Nb, V, T), the thermodynamic behaviour of the system is governed by the 
canonical partition function 

Z={{eM-mj^, (3) 

where /? = l/ksT and (• • •)^ denotes a classical canonical trace over the 
coordinates of particles of type a: 

(exp(-/?if))„ = j'^^'^' j exp(-/3i?), (4) 

with Aa being the respective thermal de Broglie wavelength. 

The two-component mixture can be formally mapped onto an equivalent 
one-component system by performing a restricted trace over the coordinates 
of only the h particles, keeping the a particles fixed. Thus, without approxi- 
mation, 

Z = (exp(-/3if„„) (exp[-/3(F66 + H,b)])b)„ = (exp(-/3ifeff))„ , (5) 

where 

HeS = Haa + Fb (6) 

is the effective Hamiltonian of the equivalent one-component system and 



-fcsTln (exp [-(}{Hbb + Hab)])^, 



(7) 
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can be physically interpreted as the Helmholtz free energy of the b particles 
in the presence of the fixed a particles. Equations (5)-(7) provide a formally 
exact basis for calculating the effective interactions. It remains, in practice, 
to explicitly determine the effective Hamiltonian by approximating the en- 
semble average in (7). Next, wc describe three general methods of attack - 
response theory, density-functional theory, and distribution function theory. 
While these methods ultimately give equivalent results, they have somewhat 
differing origins and conceptual interpretations. 

3.2 Response Theory 

The term "response theory" can have varying specific meanings, depending on 
discipline and context, but is used here to denote a collection of statistical me- 
chanical methods that describe the response, in a multicomponent condensed 
matter mixture, of the density of one component to the potential generated 
and imposed by another component. Response theory has been systematically 
developed and widely applied, over the past four decades, in the theory of sim- 
ple metals [20, 24, 25] to describe the quantum mechanical response of valence 
electron density to the electrostatic potential of metallic ions. More recently, 
similar methods have been carried over and adapted to classical soft matter 
systems, in particular, to charge-stabilised colloidal suspensions [32, 33, 34], 
polyelectrolytes [35, 36, 37], and colloid-polymer mixtures [38, 39]. Although 
most applications have been restric;tcd to the linear response approximation, 
which assumes a linear dependence between the imposed potential (cause) and 
the density response (effect), increasing attention is being devoted to nonlin- 
ear response. Below, we outline the key elements of response theory, including 
both linear and nonlinear approximations, in the context of a classical binary 
mixture. 

Perturbation Theory 

To approximate the free energy (7) of one component (6 particles) in the pres- 
ence of another component (a particles) it is often constructive to regard the 
a particles as generating an "external" potential that perturbs the b particles 
and induces their response. This external potential, which depends on the ab 
interaction and on the number density Po(r) of a particles, can be expressed 



as 





and the ab interaction term in the Hamiltonian as 




(9) 



where 
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p«(r)=^<5(r-r,) (10) 
j=i 

is the number density operator for particles of type a {a = a,b). 

If the a particles possess a property (e.g., electric charge or size) that can 
be continuously varied to tune the strength of i'ext(r)) then Ff, can be approx- 
imated via a perturbative response theory. A prerequisite for this approach is 
accurate knowledge of the free energy of a reference system of pure b particles 
(unperturbed by the a particles). Relative to this reference system, the free 
energy can be expressed as 



^0+ / 

Jo 



F, = Fo+ dX^^=Fo+ dX{H,t,)^, (11) 







J{ 
where 

Fo = FbiX = 0) = -fcBTln (exp(-/3iJ(,b))^ (12) 
is the reference free energy, 

Fb{X) = -ksTln (exp [-l3{H^b + ^H^b)])^ (13) 

is the free energy of h particles in the presence of a particles "charged" to a 
fraction A of their full strength, and 



. _ (gabexp[-/3(gbb+Agafc)]), ^ dFbjX) 

^ "''^ (exp[-/3(i/b6 + Ai?ah)])b dX ^ > 



denotes a trace of Hab over the coordinates of the b particles in this inter- 
mediate ensemble. (To simplify notation, we henceforth omit the subscript b 
from the trace over the coordinates of b particles: (• • •)^ = (• • •).) 

Applying now a standard perturbative approximation, adapted from the 
theory of simple metals [20, 24, 25], the ensemble-averaged induced density of 
b particles may be expanded in a functional Taylor scries around the reference 
system [i'ext(r) = 0] in powers of the dimensionless potential u{r) = — /3uext(r): 

{pbir)) = m + Y, /"dri ••• / dr„G("+i)(r-ri,...,r-r„)u(ri)---u(r„), 

n=i 

(15) 

where rib = ^b/V is the average density of b particles and the coefficients 

G("+^) (r - n, . . . , r - r„) = lim f -J^MllL-") (16) 



u-yQ \ du{ri) ■ ■ ■ 5u{rn] 

are the (n + l)-particlc density correlation functions [20] of the reference sys- 
tem. Equation (15) has a simple physical interpretation: the density of b par- 
ticles induced at any point r results from the cumulative response to the 
external potentials at all points {ri, . . . ,r„}, propagated through the system 
via multiparticle density correlations. 
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Further progress follows more rapidly in Fourier space, where the self- 
interaction terms in the Hamiltonian (1) can be expressed using the identity 



^l]^)aa(fc)[p„(k)p„(-k)-iV„], (17) 



while the cross-interaction term (2) takes the form 

{Hab)^=^Y.*ab{k)Pa{^){Pb{-^))^- (18) 
k 

Here Vai3{k) {a, (3 = a,b) is the Fourier transform of the pair potential Va0{r) 
and 

p«(k) = J drp«(r)e-*-- (19) 

is the Fourier transform of the number density operator (10), with inverse 
transform 

P.ir) = yJ2P-(^>'^"'- (20) 

k 

The inverse transform is expressed as a summation, rather than as an in- 
tegral, to allow the possibility of isolating the fc = component to pre- 
serve the constraint of fixed average density in the canonical ensemble: 
Pa{k = 0) = J drpa(r) = A^„- For charged systems, which interact via bare 
Coulomb pair potentials, special care must be taken to ensure that all long- 
wavelength divergences formally cancel (see Sec. 4.1 below). 
Now Fourier transforming (15), we obtain 

(P6(k)) =G(2)(/e)u(k) + -^^G(3)(k',k-k')n(k')n(k-k') + ---, kj^O, 

k' 

(21) 

where the coefficients G^"^ (Fourier transforms of G^")) are related to the 
n-particle static structure factors of the reference system via G^") = nbS^"\ 
with the static structure factors being explicitly defined by [20] 

S(2)(fc)^5(fc) = ^(p,(k)p,(-k)) (22) 

and 

5(")(ki, • • • , k„_i) = (p6(ki) • • • pb(k„_i)p,(-ki - ... - k„_i)) , n > 3. 

(23) 

Substituting M(k) = — /3?)a6(fc)yOa(k) [from (8)] into (21), the induced density 
of b particles can be expressed in the equivalent form 
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(Pb(k)) = x{k)Vab{k)Pa{^) + ^ - ^'Kbik')Vab{\k - k'|) 

k' 

X p„(k')/3„(k-k') + ---, k^O, (24) 



where 

X(fc) = -PnbS{k) (25) 
is the Unear response function and 

X'(k', k - k') = ^/3^nf,5(3) (k', k - k') (26) 

is the first nonhnear response function of the reference system. 

Substituting the equilibrium density of b particles (24) into the ab interac- 
tion (18), the latter into the free energy (11), and integrating over A, yields the 
desired free energy of the b particles to third order in the macroion density: 

Fb = Fo + NaHb lim v^b{k) + i^V] x{k) [vab{k)]' /3a(k)p„(-k) 
fc— >o zv '—^ 

+ ^ E E ^'(k'' -k - V)v^b{k)v^b{^)vabi\\^ + k'l) 

k^^O k' 

X p„(k)p„(k')/3„(-k-k'). (27) 

Finally, this free energy may be substituted back into (6) to obtain the effec- 
tive Hamiltonian. Evidently, the term in that is quadratic in yOa(k), arising 
from the term in (pb(k)) that is linear in Pa(k), is connected to an effective 
interaction between pairs of a particles. Similarly, the term in F^ that is cu- 
bic in /5o(k), coming from the quadratic term in (/3b(k)), is connected to an 
effective interaction among triplets of a particles. 



Effective Interparticle Interactions 

To explicitly demonstrate the connections between the free energy F^ and the 
effective interactions, we first identify 

4'd(fc)=x(fc)[^)a6(fc)]' (28) 

as the interaction between pairs of o particles induced by surrounding b parti- 
cles, in a linear response approximation [32, 33, 34]. As expected, the induced 
interaction depends on both the bare ab interaction and the response of the b 
particles to the external potential of the a particles. Combining the bare aa 
interaction with the induced interaction yields the linear response prediction 
for the effective pair interaction: 

=^aa(fc) + 4'dW- (29) 
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The term on the right side of (27) that is second-order in /3a(k) can be ma- 
nipulated using the identity [from (17)] 



E ^iniinj) = ^ E 4'l(fc)/5a(k)p„(-k) + iV„n„ lim - iVa^LI(O), 

(30) 

where ria = Na/V is the average density of a particles. Similarly, identifying 

t)(^)(k,k') = 2x'(k', -k - k')«a6(fc)«a6(fc')«a6(|k + k'|) (31) 

in (27) as an effective three-body interaction, arising from nonlinear response, 
and invoking the identity 



i^j^k=l k k' 

- 3p„(k)p„(-k)+27VJ, 



(32) 



the effective Hamiltonian acquires the following physically intuitive structure: 
H,s = E + -Y^ ^isinj) + ^ E ^^^(^ij'^ik), (33) 



where E, v'^g (r), and (r, r') are, respectively, a one-body "volume energy" 
and effective pair and triplet interactions, induced by the b particles, between 
the a particles. 

A natural by-product of the reduction to an effective one-component sys- 
tem, the volume energy is entirely independent of the a particle positions. 
Collecting coordinate-independent terms, the volume energy, expressed as 
E = Eiin + AE, comprises a linear response approximation. 



2 '""^ ' -fc^o 
and nonlinear corrections, the first-order correction being 



AE=^ 
6 



-fffH0,0)-^^i)f,)(k,0) . 

k 

(2) 

The effective pair interaction v^^ (r) in (33) is the transform of 



where 



^4^w = ^E^S\k,k')-fi)f,^(k,o) 



(34) 

(35) 

(36) 
(37) 
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is the first nonlinear correction to the effective pair potential, while the effec- 

tive triplet interaction w^g (r, r') is the Fourier transform of (31). Note that the 
final term on the right sides of (34), (35), and (37) arise from the constraint 
of fixed average density. 

Three observations are in order. First, since the volume energy depends, 
in general, on the mean densities of both a and b particles, it contributes 
to the total free energy and, therefore, thermodynamics of the system. This 
point has special significance in applications to charged systems, as disciisscd 
in Sec. 4.1. Second, nonlinear response of the b particles generates not only 
effective many-body interactions among the a particles, but also corrections 
to both the effective pair interaction and the volume energy. In fact, as is clear 
from (35) and (37), the nonlinear corrections to E and v^g (r) are intimately 
related to many-body interactions. Third, any influence of bb interactions on 
the effective interactions enters through the free energy and response functions 
of the reference system. Thus, the quality of the effective interactions is limited 
only by the accuracy to which the structure and thermodynamics of the pure 
6 fluid are known. 

Physical Interpretation 

While response theory is most easily formulated in Fourier space, its physical 
interpretation is perhaps more transparent in real space. The induced pair 
interaction in the linear response approximation (28) can be expressed in 
terms of real-space functions as 



Here x(|ri — r2|) is the real-space linear response function, which describes 
the change in the density of b particles induced at point r2 in response to 
an external potential applied at point ri. Referring to Fig. 1, (38) can be 
interpreted as follows. A particle of type a, centred at the origin in Fig. 1, 
generates an external potential Vab{ri), which acts on b particles at all points 
ri. This potential induces at point r2 a change in the density of b particles 
given by J dri x(|ri — r2\)vab{ri). This induced density, which depends on pair 
correlations (via x) i^i the intervening b fluid, then interacts with a second 
a particle, at displacement r from the flrst. The net result is an effective 
interaction between the pair of a particles induced by the medium. The linear 
response contribution to the volume energy (per particle) associated with ab 
interactions (34) has a closely related form: 



The physical interpretation is similar, except that the induced density now 
interacts back with the flrst a particle, generating a one-body (self) energy. 
An analogous interpretation applies to nonlinear response and induced many- 
body interactions [40]. 




dri / dr2x(|ri - r2|)i;a6(ri)z;ab(|r2 - r|). 



(38) 




(39) 
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3.3 Density-Functional Theory 

An alternative, yet ultimately equivalent, approach to deriving effective in- 
teractions, follows from classical density- functional theory (DFT) [41, 42]. 
Classical DFT has a long history, dating back a half-ccntury to the earliest 
integral-equation theories of simple liquids [20]. Following the establishment 
of formal foundations [43] , DFT has been widely applied, in recent decades, to 
a variety of soft condensed matter systems, including colloids, polymers, and 
liquid crystals. Connections between density-functional theory and effective 
interactions in charged colloids have been established by Lowen et al [44, 45] 
and van Roij et al [46, 47]. The essence of the theory is most easily grasped 
in the context of an ab mixture, where the challenge again is to approximate 
the free energy (7) of a fluid of b particles in the presence of fixed a particles. 

The basis of the density-functional approach is the existence [43] of a 
grand potential functional i7h[p(,] the square brackets denoting a functional 
dependence - with two essential properties: f2b[Pb] is uniquely determined by 
the spatially- varying density Pbi'^), for any given external potential Woxt(r), 
and is a minimum, at equilibrium, with respect to Pb(r). The grand potential 
functional is related to the Helmholtz free energy functional Tb[Pb\ via the 
Legendre transform relation 

^b[pb] = :Fb[pb] - HbNb, (40) 

where fXb is the chemical potential of b particles. The free energy functional 
naturally separates, according to 

^b [Pb] = ^id [Pb] + ^ext [Pb] + ^ex [Pb] , (41) 

into an ideal-gas term .Fid, which is the free energy in the absence of any 

interactions, an "external" term J-cxt, which results from interactions with 
the external potential, and an "excess" term J^ex, due entirely to interparticle 
interactions. The purely entropic ideal-gas free energy is given exactly by 

^id[P6] = ksT J drp6(r)[ln(pb(r)4) - 1], (42) 

where Ab is the thermal wavelength of the b particles, while the external free 
energy can be expressed as 

•^ext[P6] = j drp6(r)vext(r), (43) 

which is equivalent to Hab in (9). Inserting (42) and (43) into (41) yields 

n[pb] = kBT j drpfe(r)[ln(/)b(r)ylf)-l] + J dr pbir)v^^t{r) + J'eM (44) 



The excess free energy can be expressed in the formally exact form. 
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^e4Pb] j j J ^^P'ni^'"bb;r,r']vbb{\r-r'\) 

= 11^"^^ J dv'pb{v)pb{v') d\gf^[Xvbb;v,v']vbb{\v-v'\), (45) 

where pj^^^ [Ai'^^; r, r'] is the two-particle number density (a unique funetional 

('21 

of the pair potential), g)^^'[\%Jbb'-,^-,^'] if' the corresponding pair distribution 
functional, and A is a couphng (or charging) constant that "turns on" the 
intcrparticlc correlations. In general, the pair distribution functional is not 
known exactly and must be approximated. For weakly correlated systems, it 
is often reasonable to adopt the mean-field approximation 

= \j '^'^ j <^^' Pb{r)pb{r')vbb{\v - r'l), (46) 

which amounts to entirely neglecting correlations and assuming gl^^' = 1. 

In a further approximation, valid for weakly inhomogeneous densities, 
the ideal-gas free energy functional is expanded in a functional Taylor series 
around the average density, which is truncated at quadratic order: 

/idN ^ fid(»i) + (j^) / [^'(f) - "'1 

The first term on the right, Fid(n6) = NbkBT^n{nbAb) — l], is the ideal-gas free 
energy of a uniform fluid of h particles. The second (linear) term on the right 
vanishes identically by virtue of the constraint of constant average density. 
Evaluating the functional derivative in the of the third (quadratic) term, 

^^^idN _ fci^^.5(r_r'), (48) 



5p{r)5p{r') pb{r) 

and combining (41), (43), (46), and (47), the mean-field free energy functional 
finally can be approximated by 

J'bipb] - NbkBT[ln{nbAl) - 1] + ^ dr [/3fc(r) - Ub]^ + J dr pb{r)veMr) 

+ Av' pb{v)pb{v')vbb{\v - r'l). (49) 

The equilibrium density is now determined by the minimization condition 

p5Qt\p^ = ln(n(,4)+^-l+/3i;ext(r)+/3 / dv' pb{v')vbb{\r-v'\)-(3nb = 0. 
opb{r) ^ rib J 

(50) 
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Fourier transforming and solving for the equilibrium density yields 
(k)> = = X(fc)^)ext(k), k ^ 0, 



where 



+ l3nbVbb{k) 



X{k) 



-j3nb 



(51) 



(52) 



1 + 0nbVbb{k) 

is a mean-field approximation to the linear response function introduced above 
in (25). Now expressing the free energy functional (49) in terms of Fourier 
components, 



J^b[pb] ^ NbkBT[ln{nbAl) -l]+nb lim 



fc-»0 



^NbVbbik) + NaVabik) 



+ V I]^f(k)^a6(fc)Pa(k) + ^12 (*'>''(^) + Pb(k)pf,(-k), (53) 

and substituting for the equilibrium density from (51), we obtain - to second 
order in the a particle density - the equilibrium Helmholtz free energy of b 
particles in the presence of the fixed a particles: 



Fb = NbkBT[ln{nbAl) - 1] + rib Hm 



fc-»0 



^NbVbbik) + NaVabik) 



(54) 



After identifying 



Fo = NbkBT[ln{nbAl) - 1] + -Nbm lim Vbb{k) 

A fe— »0 



(55) 



as the free energy of the uniform h fluid in the absence of a particles, (54) 
is seen to have exactly the same form as (27), to quadratic order in /5a(k). 
The same effective interactions thus result from linearized classical density- 
functional theory as from linear response theory. Moreover, the same agree- 
ment is also found for the effective triplet interactions derived from nonlinear 
response theory [40] and nonlinear DFT [48]. 



3.4 Distribution Function Theory 

Still another statistical mechanical approach to calculating eff'ective interac- 
tions is based on approximating equilibrium distribution functions. This ap- 
proach has been developed by many workers and applied to charged colloids 
in the forms of various integral-equation theories [49, 50, 51, 52, 53, 54, 55] 
and an extended Dcbyc-Hiickcl theory [56, 57, 58, 59]. The basic elements of 
the method are sketched below, again in the context of a simple ah mixture. 
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As shown above in Sec. 3.1, the key quantity in any theory of effective 

interactions is the trace over degrees of freedom of the b particles of the ab 
interaction term in the Hamiltonian [see (2) and (11)]. This partial trace can 
be expressed in the general form 



{Hab) = J^^J dr' (Pa(r)p6(r'))^;a6(|r-r'|), 



(56) 



where (• • •) denotes an ensemble average over the coordinates of the b particles 

with the a particles fixed. The density of the fixed a particles being unaffected 
by the partial trace, we can replace {pa{r)pb{r')) by Pa(r) (pb(r')) in (56): 

{Hab) = "f-Jdr 1 dr'pa(r)5i6^(r - r')t;a6(|r - r'|), (57) 

(2) 

thus introducing the ab pair distribution function (r), which is defined 

(2) 

via (p6(r)) ~ nbg)j^fl {v). This pair distribution function is proportional to the 
probability of finding a b particle at displacement r from a central a particle. 

More precisely, given an a particle at the origin, ribg^^ (r)dr represents the 
average number of b particles in a volume dr at displacement r. 

Distribution function theory evidently shifts the challenge to determining 
the cross-species [ah] pair distribution function. To this end, wc consider an 
approximation scheme rooted in the theory of simple liquids [20] that 
illustrates connections to response theory and density-functional theory. The 
starting point is the fundamental relation, valid for any nonuniform fluid, 
between the equilibrium density, an "external" applied potential, and the 
(internal) direct correlation functions. Minimization of the Helmholtz free 
energy functional (41) with respect to the density, at fixed average density, 
yields the Euler-Lagrange relation [43] 

<Pfe(r)) = -73- exp(-/3i)ext(r) + ^ [Pb\ rj), (58) 

where the one-particle direct correlation functional (DCF), defined as 

c«[p6;r] = -fi^-j^ = -/3Mex[p6;r], (59) 

is a \miquc functional of the density that is proportional to the excess chemical 
potential of b particles /iex (associated with bb interparticle interactions). 

Although it provides an exact implicit relation for the equilibrium density, 
(58) can be solved, in practice, only by approximating c^^'' [p;,; r]. Approxima- 
tions are facilitated by expanding c'^^ [pb] r] in a functional Taylor series about 
the average (bulk) density Ub'- 

c\^\pb-A=c^i^\nb)+ fdr'cg\r-r';nb)[{pb{r'))-nb] + ---, (60) 
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where 

4'(r-r;n,)=^^hm^^(^^-^j (61) 

is the two-particle DCF of the (reference) uniform fluid and, more generally, 

(") r 1 (5""^4^^[/'b;i"i] u-n\ 

ct-.j, Pb;ri,...,r„ = . . . — . , X 62 
<5p6(r2) • • •(5pb(r„) 

is the n-particle DCF. Note that higher-order terms in the series, correspond- 
ing to multiparticle correlations, are nonlinear in the density. 

Substituting (60) into (58), and identifying (e^'"'/^g) exp[4^^(n6)] as the 
bulk density rib, the nonuniform equilibrium density can be expressed as 

(P6(r)) = n^exp ^-/J'!;ext(r) j 6.t' cf^ {t - r' ■,nb)[{pb{Y')) - nb] + ■ ■ 

(63) 

Various practical approximations for the DCFs correspond, within the frame- 
work of integral-equation theory, to different closures of the Ornstein-Zernike 
relation [20], 



h''^{v)=c''^{v)+nb j dr'cf,)(r-r')/if,)(r'), 



(64) 



which is an integral equation relating the two-particle DCF to the pair corre- 
lation function, /i^^''(r) = ''(r) — 1- 

Truncating the functional expansion in (63) and retaining only the term 
linear in density, thus neglecting multiparticle correlations in a mean-field 
approximation, is equivalent to the hypernetted-chain (HNC) approximation 
in integral-equation theory. Making a second mean-field approximation by 
equating the two-particle DCF c\^'{T;nb) to its asymptotic limit [20], 

lim cf^{Y;nb) = -Pvbb{r), (65) 

r— >oo 

thereby neglecting short-range correlations, (63) reduces to 

(P6(r)) = ntexp [-I3v,^{v) - (3 J dr' vil\\r - r'\)[{pb{r')) - nb]j , (66) 

corresponding to the mean- spherical approximation (MSA) in integral-equation 

theory. In passing, we note that (66) also provides a basis for the Poisson- 
Boltzmann theory of charged colloids and polyelectrolytes, if we identify 

96V(r) = ^;ext(r) + J dr' ^;(^)(|r - r'\)[{pb{r')) - Ub] (67) 

as the mean-field electrostatic potential energy of a 6 particle (charge qb) 
brought from infinity, where the potential tp ~ 0, to a displacement r away 
from an a particle, and combine this with Poisson's equation, 
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rib] ■ (69) 



V^^Piy) = --^Y.1iP,(y), (68) 

i 

where e is the dielectric constant of the medium and the sum is over all species 
with charges qi and number densities pi{Y). 

If we now make one further approximation by linearizing the exponential 
function, valid in the case of potential energies much lower than thermal 
energies, then (66) becomes 

(p,(r)) ^n,[l- pv^r) - [3 j dv' v^^ {\v - v' \)[{p,{v' )) 

Fourier transforming (69) and solving for the equilibrium density finally yields 

1 + pnbVbb[K) 

which is identical in form to the linear response and DFT predictions [sec (24) 
and (51)] with the same linear response function x(/e) as before [see (52)]. 

Alternatively, we may first linearize the exponential in (66) and then ex- 
ploit the Ornstcin-Zcrnike relation (64) to solve recursively for the equilibrium 
density, with the result 

{Pb{v)) =nb(l- /3t;ext(r) + J dr'4^^(r - r')[-(3nbV,^t{r') + ■ ■ 

= nb- pUb J dr'[(5(r') + Ubh^^^ {r')]v,^,{r'), (71) 

where 

-(3v4Sir)+nbh^\r)]=xir) (72) 

can be identified as the real-space linear response function. By Fourier trans- 
forming (71) and (72), we recover the linear response relation (70) with a 
linear response function, 

X(fc) = -f3nb[l + nbhf^j(k)] = -pHbSik), (73) 

defined precisely as originally in (25). The connection to the linear response 
function in (52) is established via the Fourier transform of the Ornstein- 
Zernike relation (64), 

= . -^^g'g' . ,74, 

l-ntcg'(fc) 1 + /Jn,«g'(l!) 

where we have assumed a mean-field (random phase) approximation for the 

two-particle DCF, c!j^{k) ^ —Pi'bb'i^-) ["^^ (65)]- Distribution function theory 
therefore predicts the same linear response of the 6-particle density to the 
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a-particle external potential - and thus the same effective interactions - as do 
both response theory and density-functional theory. 

In closing this section, we emphasize that the effective interactions derived 
here in the canonical ensemble apply to experimental situations in which the 
particle densities arc fixed. In many experiments, however, the system may be 
in chemical equilibrium with a reservoir of particles, allowing fluctuations in 
particle densities. The appropriate ensemble then would be the semigrand or 
grand canonical ensemble for a reservoir containing, respectively, one or both 
species. It is left as an exercise to derive the effective interactions in these 
other ensembles. (Hint: Consider carefully the appropriate reference system.) 

4 Applications 

The preceding section describes several general, albeit formal, approaches to 
modelling effective interparticle interactions in soft matter systems. To illus- 
trate the practical utility of some of the methods, the following section briefly 
outlines applications to three broad classes of system: (1) charge-stabilised col- 
loidal suspensions, (2) colloid-polymer mixtures, and (3) polymer solutions. 
These systems exhibit two of the most common forms of microscopic interac- 
tions, namely, electrostatic and excluded-volume interactions. Further details 
can be found in [28, 33, 34, 38, 39, 40]. 

4.1 Charged Colloids 
Primitive Model 

Charge-stabilised colloidal suspensions [7, 8, 9, 10] are multicomponent mix- 
tures of macroions, counterions, and salt ions dispersed in a molecular solvent 
and stabilised against coagulation by electrostatic interparticle interactions. 
For concreteness, we assume aqueous suspensions (water solvent). A reason- 
able model for these complex systems is a collection of charged hard spheres 
and point microions interacting via bare Coulomb pair potentials in a dielec- 
tric medium (Fig. 2). The point-microion approximation is valid for systems, 
such as colloidal suspensions, with large size asymmetries between macroions 
and microions. The macroions, of radius a (diameter a = 2a), are assumed to 
carry a flxed, uniformly distributed, surface charge —Ze (valence —Z), which 
may be physically interpreted as an effective charge, renormalized by associ- 
ation of oppositely charged counterions (valence z <^ Z).ln & closed system, 
described by the canonical ensemble, global charge neutrality constrains the 
macroion and counterion numbers, and A^c, via the relation ZN^ = zNc- 

The salt-water mixture is modelled as an electrolyte solution of Ng dis- 
sociated pairs of point ions of valences ±z. The microions then number 
A''+ = iVc -I- Ns positive and A_ = negative, totaling = Nc + 2Ns. 
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Fig. 2. Left: Primitive model of a charge-stabilised colloidal suspension, consisting 
of hard spherical macroions (valence diameter a) and point monovalent microions 
(counterions, salt ions) dispersed in a dielectric continuum. Right: Effective one- 
component model of pseudomacroions, after coarse-graining of microions. 

Within the coarse-grained "primitive" model of charged colloids, the wa- 
ter is treated as a dielectric continuum, characterized entirely by a dielec- 
tric constant e. This approximation amounts to preaveraging over the vast 
number of solvent degrees of freedom. For simplicity, we completely neglect 
charge-induced-dipole and other polarization interactions [60, 61, 62], which 
are shorter-ranged than charge-charge interactions and vanish if solvent and 
macroions are index-matched (i.e., have the same dielectric constant). The 
bare electrostatic interactions can bo represented by Coulomb pair potentials: 
Vmm{r) = Z'^f? jtr (r > cr), Vccir) = z^e^/er, and Vmc{r) = Zze^/er (r > a). 
Note that in the primitive model, the solvent acts only to reduce the strength 
of Coulomb interactions by a factor 1/e. In addition, the macroion hard cores 
interact via a hard-sphere pair potential. 

Response Theory for Electrostatic Interactions 

Following the methods of response theory laid out in Sec. 3.2, effective interac- 
tions now can be derived by integrating out the degrees of freedom of the mi- 
croions, reducing the multicomponent mixture to an effective one-component 
system of pseudo-macroions [32]. To simplify the derivation, we first con- 
sider the rather idealized case of salt-free suspensions. The bare Hamilto- 
nian of the two-component model then decomposes naturally, according to 
H = Hmm + Hcc + Hmc, into & macroion term 



where iJns is the Hamiltonian for neutral hard spheres (macroion hard cores), 
a counterion term 




(75) 




(76) 
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where K^. is the counterion kinetic energy, and a macroion-counterion inter- 
action term 

Nm Nc 

Hmc = XI ^'^a ) • (^^) 

1=1 j=l 

By analogy with (11), the free energy of the counterions in the external 
potential of the macroions can be expressed as [20, 32]: 



Fr. = Fn 



mc/x > 



(78) 



where Fq = — fcsTln (exp(— /3iJc)) is now the reference free energy of the coun- 
terions in the presence of neutral (hard-core) macroions, and the A-integral 
adiabatically charges the macroions from neutral to fully charged. Neglecting 
counterion structure induced by the macroion hard cores, neutral macroions 
would be surrounded by a uniform "sea" of counterions. As the macroion 
charge is turned on, the counterions respond, redistributing themselves to 
form a double layer (surface charge plus neighbouring counterions) around 
each macroion. 

It is a special property of Coulomb-potential systems that, because volume 
integrals over long-ranged 1/r potentials diverge, each term on the right side 
of (78) is actually infinite. Although the infinities formally cancel, it proves 
convenient still to convert Fq to the free energy of a classical one-component 
plasma (OOP) by adding and subtracting the (infinite) energy of a uniform 
compensating negative background 



^bg = —NcTlc hm Vcc{0) 
Z fe— ►o 



(79) 



where ric is the average density of counterions in the volume unoccupied by 
the macroion cores. Because the counterions are strictly excluded (with the 
background) from the hard macroion cores, the OCP has average density 
Tic = ^c/[y{^ — v)]j where rj = ^{Nm/V)a^ is the macroion volume fraction 
and V{1 — T]) is the free volume. Thus, 



Fc = Focp + / dA {Hmc) 



Jo 



bgi 



(80) 



where -Focp = -Fo + E[,g is the free energy of the "Swiss cheese" OCP in the 
presence of neutral, but volume-excluding, hard spheres. 

All of the formal expressions derived in Sec. 3.2 for a generic two- 
component [ah) mixture now carry over directly, with the identifications 
a ^ m and 6 <-> c. In the linear response approximation [33, 34], the vol- 
ume energy is given by 



^lin = Focp + ^^^L'd(O) + ^m^c lim 



■ik) 



(81) 
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the effective (electrostatic) pair potential by 



with induced potential 



(82) 



(83) 



and the effective triplet potential by 

(k, k') = 2x'(k', -k - k')Vmc{k)Vmc{k')Vmc{\^ + k'|), (84) 

where x{k) and x'{k) are linear and first-order nonlinear response functions of 
the uniform OCP. Similarly, the first-order corrections for nonlinear response 
are given by 



AE 



5]0fff)(k,k')-7V„^«^^^(k,0) 



k,k' 



(85) 



and 



^*e«^(^) = ^E<ff(k>k')-^i)S'(k,0), 



(86) 



where Vf = V{1 — rj) is the free volume. 



Random Pheise Approximation 

Further progress towards practical expressions for effective interactions re- 
quires specifying the OCP response functions. For charged colloids, the OCP 
is typically weakly correlated, characterized by relatively small coupling pa- 
rameters: r ~ Xs/ac 1, where \b = flz'^e'^/e is the Bjcrrum length and 
ttc = {S/AnncY^^ is the counterion-sphere radius. For example, for macroions 
of valence Z = 500, volume fraction r] = 0.01, and monovalent counteri- 
ons suspended in salt-free water at room temperature (Ag = 0.714 nm), 
we find F ~ 0.02. For such weakly-correlated plasmas, it is reasonable - 
at least as regards long-range interactions - to neglect short-range corre- 
lations. We can thus adopt a random phase approximation (RPA), which 
equates the two-particle direct correlation function to its exact asymptotic 
limit: c^'^\r) = —/3vccir) or S-'^\k) = — 47r/3^;^e^/efc^. Furthermore, we ig- 
nore the influence of the macroion hard cores on the OCP response functions, 
which is reasonable for sufficiently dilute suspensions. Within the RPA, the 
OCP (two-particle) static structure factor and linear response function take 
the analytical forms 
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and 

X(fc) = -/3ne5(fc) = ^^g^, (88) 



where k = ^imcz'^e^ / tkuT is the Debye screening constant (inverse screen- 
ing length), which governs the form of the counterion density profile and of the 
screened effective interactions. In the absence of salt, the countcrions are the 
only screening ions. The macroions themselves, being singled out as sources 
of the external potential for the counterions, do not contribute to the density 
of screening ions. Fourier transforming (88), the real-space linear response 
function takes the form 

X(r) = -/3ne [5{v) + n,K,{r)\ , (89) 

where 

hco{r) = -^^^ (90) 
e r 

is the coimtcrion-countcrion pair correlation function. Note the scrcened- 
Coulomb (Yukawa) form of hcc{r), with exponential screening length k^^. 
Equation (89) makes clear that there are two physically distinct types of coun- 
terion response: local response, associated with counterion self correlations, 
and nonlocal response, associated with counterion pair correlations. 

Proceeding to nonlinear response, we first note that the three-particle 
structure factor obeys the identity 

S(3)(k,k') = 5(fc)5(fc')5(|k + k'|) \l + nl c(3Hk,k')l , (91) 

where c'^^^(k, k') is the Fourier transform of the three-particle DCF. Within 
the RPA, however, c'^^^ and all higher-order DCF's vanish. Thus, from (25), 
(26), and (91), the first nonlinear response function can be expressed in Fourier 

space as 

X'(k, k') = -'^^{k)x{k')x{\k + k'l) (92) 
and in real space as 

X'(ri -r2,ri -rs) = -^^1^ J drx(|ri - r|)x(|r2 - r|)x(|r3 - r|). (93) 



Counterion Density Profile 

An explicit expression for the ensemble-averaged counterion density is ob- 
tained by substituting the RPA linear response function (88) into the linear 
response relation 

(/},(k)) = x(fc)«cxt(fc) = x{k)vmc{k)PmO<-), k ^ 0. (94) 



Inverse transforming (94) yields 
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Pe(r)=X;Po(|r-Ri|)= [dr'xi\r-r'\)Y,Vmc{\r'-Ri\), (95) 

which is the real-space linear response counterion density in the presence of 
macroions fixed at positions Ri, expressed as a sum of single- macroion coun- 
terion density orbitals po{r) - the inverse transform of po{k) = x(^)^mc(^)- 
Substituting (89) and (90) for the real-space RPA linear response function 
into (95), the linear response counterion density profile can be expressed as 



-Vmc{\r - Ril) + ^ / dr' J^Vmci\r' - Ri|) 

47r J |r — r I 



(96) 



where the two terms on the right correspond again to local and nonlocal 
counterion response. 

For hard-core macroions, the form of the macroion-coTintcrion interaction 
inside the core is arbitrary and can be specified so as to minimize counterion 
penetration inside the cores [46]. Thus, assuming 

r -Zze^ ^ > „ 



Vmc{r) = < ^ (97) 

t -^§^a, r<a 

leaves the freedom to choose the parameter a appropriately. As shown in [33] 
and [46] , at the level of linear response, penetration of counterions inside the 
macroion cores is eliminated by choosing a = Ka/{1 + ko). This choice yields 

,,, 47rZ2e^ r ,, , . /, n1 

Vmcik) = — — r-r^ cos(A;a) -I- - sm(A:a) (98) 

e(l -h K,a)k^ L k \ 

and 

z 3^ l + na r ' ^ > (99) 
0, r < a, 

which agrees precisely with the asymptotic (r ^ oo) expression predicted by 
the DLVO theory of charged colloids [22, 23]. 

Effective Electrostatic Interactions 

Practical expressions for the effective electrostatic interactions arc obtained 
by explicitly evaluating inverse Fourier transforms. By combining (81), (83), 
(84), (85), and (98), the volume energy can be expressed as the sum of the 
linear response approximation, 

E.. = Foc.-nJ-^-^-^, (100) 
2e 1 -h Aca 2 

and the first nonlinear correction. 
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= -^^^^ (/ dr [p,{rf - J dr [po{r)f^ . (101) 

The first term on the right side of (100) accounts for the counterion entropy 
and the second term for the macroion-counterion electrostatic interaction en- 
ergy. The latter term happens to be identical to the energy that would result if 
each macroion's countcrions were all concentrated at a distance of one screen- 
ing length {k~^) away from the macroion surface. 

Prom (82), (83), (88) and (98), the linear response prediction for the ef- 
fective pair interaction is given by 

which is identical to the familiar DLVO screened-Coulomb potential in the 
dilute limit of widely separated macroions [22, 23], while (86) yields the first 
nonlinear correction 

M'^ir) = J dr'po(r')po(|r - r'|) [po{\r - r'|) - ^] . (103) 

Finally, from (84) and (92), the eff'ective triplet interaction is 



(3) 



f'(ri2,ri3) = -^ /drpo(|ri-r|)po(|r2-r|)po(|r3-r|). (104) 



Note that the final terms on the right sides of (100), (101), and (103) originate 
from the charge neutrality constraint. 

The above results generalize straightforwardly to nonzero salt concentra- 
tion. Here we merely sketch the steps leading to the final expressions, referring 
the reader to [34] and [40] for details. Assuming fixed average number density 
(in the free volume) of salt ion pairs, = Ng/Vf, the total average microion 
density is = n++n_ = nc+2Ti.,, where n± are the average number densities 
of positive/negative microions. Following [34], the Hamiltonian generalizes to 
H = Hmm + + Hm+ + Hm-, whcrc is the Hamiltonian of all niic;roions 
(counterions and salt ions) and Hm± are the electrostatic interaction energies 
between macroions and positive/negative microions. The perturbation theory 
proceeds as before, except that the reference system is now a two-component 
plasma. The presence of positive and negative microion species entails a pro- 
liferation of response functions, Xij and x'ijk^ ^'J' k — ±, and a generalization 
of (89) to 

X++{r) = -I3n+ [S{r) + n+h++{r)] 

X+-(r) = -/3n+n_/i+_(r) (105) 

X— (r) = -/3n_ [^(r) + n_/i__(r)] , 

where hij (r ) , i, j = ± , are the two-particle pair correlation functions of the mi- 
croion plasma. Substituting the ensemble-averaged microion number densities 
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into the multi-component Hamiltonian yields expressions for the macroion- 

microion interaction terms and, in turn, the effective interactions. 

The ultimate effect of salt is to modify the previous results as follows. 
First, the average counterion density that appears in the Debye screening 

constant and in the linear response function (88) is replaced by the total 



average microion density: k = -J Aimfj,z'^e^ / eksT and x(^) = —(}n^S(k). The 
first nonlinear response function retains its original form (92), but with the 
new definition of k and with Uc replaced by n^. Second, the linear response 
volume energy becomes [34] 

F -F N ^^'^ ^ {N+-N.f kBT 

^-lin - ^plasma l^m I + N+ + N_ 2 ' 



where 



i^piasma = kBT{N+[\n{n+AX) - 1] + N_[Hn_At) - 1]} (107) 



is the free energy of the unperturbed microion plasma (in the free volume) 
and A± denotes the thermal wavelengths of positive and negative microions. 
Third, the effective triplet interaction and nonlinear corrections to the effective 
pair interaction and volume energy are generalized as follows: 

PAE = -^ ^^^"3"-^ (/ dr [poirt -n.jdv [p,{r)f^ (108) 



PAv^^{r) = -y^±^^ y drVo(r')po(|r-r'|)[po(|r-r'|)-!|i (109) 



/3^^ff^(ri2,ri3) = - ^"+ 3"'^ /drpo(|ri-r|)po(|r2-r|)po(|r3-r|). (110) 

These results imply that nonlinear effects increase in strength with increasing 
charge and concentration of macroions and with decreasing salt concentration, 
and that effective triplet interactions are consistently attractive. It is also clear 
that in the limit of zero macroion concentration (rig = n+ — n_ ^ 0), or of 
high salt concentration (n^ — > 00), such that (n+ — n-)/nn — > 0, the leading- 
order nonlinear corrections all vanish. This result - a consequence of charge 
neutrality - may partially explain the remarkably broad range of validity of 
DLVO theory for suspensions at high ionic strength. 

The wide tunability of the effective electrostatic interactions leads to rich 
phase behaviour in charge-stabilised colloidal suspensions. Simulation stud- 
ies [63] of one-component systems interacting via the screened- Coulomb pair 
potential have demonstrated that variation in the Debye screening constant - 
corresponding in experiments to varying salt concentration - can acc;ount for 
the observed cross-over in relative stability between stable fee and bee crystal 
structures. The density dependence of the volume energy, resulting from the 
constraints of fixed density and charge neutrality, can have profound implica- 
tions for thermodynamic properties (e.g., phase behaviour, osmotic pressure) 



Effective Interactions in Soft Materials 



27 



of highly deionized suspensions. Specifically, the volume energy has been pre- 
dicted by DFT [46, 47], extended Dcbyc-Hiickcl theory [59], and response 
theory [64] to drive an unusual counterion-induced phase separation between 
macroion-rich and macroion-poor phases at low (sub-mM) salt concentrations. 
Despite purely repulsive pair interactions, which oppose bulk phase separa- 
tion, counterion entropy and macroion self-energy can, according to predic- 
tions, conspire to drive a spinodal instability. It remains unresolved, however, 
whether this predicted instability is related to experimental observations of 
anomalous phase behaviour in charged colloids, including stable voids [65] and 
metastable crystallites [66]. 

Predictions of response theory can be directly tested against simulations. 
As an example. Fig. 3 presents a comparison [40] with available data from 
ab initio simulations [67] for the total potential energy of interaction between 
a pair of macroions, of diameter a = 106 nm and valence Z = 200, in a 
cubic box of length 530 nm (with periodic boundary conditions) in the ab- 
sence of salt. The theory is in excellent agreement with simulation, although 
nonlinear effects are relatively weak for these parameters. Figure 3 also illus- 
trates the effective triplet interaction between a trio of macroions arranged in 
an equilateral triangle for <t = 100 nm and two different valences, Z = 500 
and 700, computed from (110) [40]. The strength of the attractive interaction 
grows rapidly with increasing macroion valence and with decreasing separation 
between macroion cores. Other methods, including DFT [48] and Poisson- 
Boltzmann theory [68], predict qualitatively similar triplet interactions. In 
concentrated suspensions of highly- charged macroions, higher-order effective 
interactions may become significant. 

4.2 Systems of Hard Particles 

Soft matter systems that contain hard (impenetrable) particles include col- 
loidal and nanoparticle dispersions and coUoid-nanoparticle mixtures. In such 
systems, the bare interparticle interactions depend, at least in part, on the 
geometric volume excluded by each particle to all other particles. In general, 
van der Waals, electrostatic, and other interactions may also be present. The 
simplest case is that of spherical particles, as in bidisperse or polydisperse mix- 
tures of colloids and/or nanoparticles. In a binary {ab) hard-sphere mixture, 
the bare pair interactions have the form 

r \ f oo, r < Ra + Rb riii\ 
^«'^W = |o, r>R^ + R;, (111) 

where a,f3 = a,b and Ra denotes the radius of particles of type a. For other 
shapes, the pair interactions naturally depend also on the orientations of the 
particles. The methods outlined in Sec. 3 for modelling effective interactions 
are easily adapted to hard-particle systems. In what follows, we first present 
general formulae and then describe an application to a common model of 
colloid-polymer mixtures. 
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Fig. 3. Left: Interaction energy of two macroions (diameter 106 nm, valence 200) 

in a cubic box of length 530 nm (with periodic boundary conditions) at zero salt 
concentration. The potentials are shifted to zero at maximum macroion separation. 
Dashed curve: linear response prediction. Solid curve: nonlinear response predic- 
tion [40]. Symbols: ab initio simulation data [67]. Right: Effective triplet interac- 
tion between three macroions, arranged in an equilateral triangle of side length r, 
with macroion diameter a = 100 nm, macroion valence Z = 500 (dashed curve) or 
Z = 700 (solid curve), volume fraction = 0.01, and salt concentration Cs = 1 jjM. 
Computed from (110) [40]. 

Response Theory for Excluded- Volume Interactions 

As noted in Sec. 3.2, the perturbative response theory is most useful in cases 
where the mesoscopic particles possess a property that can be continuously 
varied to turn on the "external" potential for the other (microscopic) particles. 
In the case of electrostatic interactions, the obvious timable property is the 
charge on the particles. Analogously, for excluded-volume interactions, the 
relevant property is the volume occupied by the particles. In an ab mixture, 
we can imagine the mesoscopic (a) particles to be inflated continuously from 
points to their full size. To grow in size, the a particles must push against the 
surrounding b particles. As they grow, the a particles sweep out spheres, of 
radius Ra + Rt, from which the centres of the b particles arc excluded. The 
change in free energy during this process equals the reversible work performed 
by the expanding a particles against the osmotic pressure exerted by the b 
fluid. Adapting (11) to quantify this conceptual image, the Helmholtz free 
energy of the b particles in the presence of the a particles can be expressed as 



where Fq is again the free energy of the unperturbed (reference) b fluid, the 
A integral continuously scales the a particles from points to full size, and 
KxcIPalA] and Ublpa, Pbl X] - both functionals of density - are, respectively, 
the total volume excluded to the b particles, and the osmotic pressure exerted 
by the b fluid, in the presence of a particles expanded to a fraction A of their 
full volume. 




(112) 
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Although (112) is formally exact for any shape of particle, the compli- 
cated dependence of the excluded volume and osmotic pressure on the scale 
parameter A precludes an exact evaluation of for an arbitrary configuration 
of a particles. One approximation scheme is based on expanding the osmotic 
pressure in powers of A. Expanding around A = 0, for example, yields 

-^j dAAKxc[/Oa;A] + ---, (113) 

where nf'\nh) is the osmotic pressure of the reference fluid and Kxc[Pa] is the 
total volume excluded to the h particles by full-sized a particles. In practical 
applications, such as colloid-polymer mixtures, the system is often coupled 
to an infinite reservoir of b particles (e.g., a polymer solution) that fixes the 
chemical potential (ib of the h component. In this case, a more natural choice 
for il^"' may be the osmotic pressure of the reservoir leading to the 

approximation 

F6 = Fo + 77^(''Vexc[Pa]. (114) 

The higher-order terms in (113), which depend on hh pair interactions, are dif- 
ficult to evaluate and commonly ignored. Nevertheless, (114) is often a reason- 
able approximation, especially in the dilute limit (ria — >■ 0), where iTf,[/9o, p^] A] 
depends weakly on Ua and A. In the idealized case in which the h fluid can 
be modelled as a noninteracting gas, the osmotic pressures of the system and 
reservoir are equal and (114) then becomes exact (see below). 

It remains a highly nontrivial problem, for an arbitrary conflguration of a 
particles, to approximate the excluded volume KxciPa], which requires deter- 
mining the intersection volumes of the mutually overlapping exclusion spheres 
that surround the a particles. In the simplest case of non-overlapping spheres, 
the total excluded volume is merely the excluded volume of a single a particle 
times the number of particles. In the next simplest case, in which only pairs 
of exclusion spheres overlap, the intersection volumes of overlapping pairs 
must be subtracted to avoid double-counting. The next step corrects for the 
intersection volume of three mutually overlapping spheres. The geometrical 
problem is illustrated in Fig. 4 for a mixture of hard spherical colloids and 
coarse-grained spherical polymers. Expressing the total excluded volume as a 
sum of overlap terms, and assuming isotropic interactions (spherical particles), 
the free energy can be approximated by 

\ 1^3 = 1 ■ i^j^k=l J 

(115) 

where Ve^c is the excluded volume of a single a particle, Vov\rij) is the 
intersection volume of a pair of overlapping exclusion spheres surrounding 
particles i and j, and Vm^ (r^ , r^^) is the intersection volume of three mutually 
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overlapping spheres surrounding particles i, j, and k. In the more general case 
of anisotropic interactions (nonspherical particles), the effective interactions 
depend also on the relative orientations of the particles. 




Fig. 4. Left: Colloid-polymer mixture with hard (excluded-volume) interactions. 

Depletion of polymer coils from spaces between colloids induces effective interactions 
between colloids. Right: Asakura-Oosawa-Vrij model with colloids treated as hard 
spheres (diameter tJc) and polymers as coarse-grained spheres (diameter ap). The 
effective interactions are a function of the total volume excluded by the colloids to the 
polymer centres, which depends on the intersection volumes of mutually overlapping 
spheres of exclusion (diameter ctc + Op) surrounding each colloid. 

Prom (115), the effective Hamiltonian is seen to have the same general 
form as in (33), with a one-body volume term 

an effective pair potential 

^^^^W==«.a(r)-<V„(,2)(r), 
and an effective triplet potential 

v^^i^ij^^ik) = ni'-^V(^\v,j,Vik). (118) 

This simple low-density approximation exhibits several noteworthy features. 
The induced pair potential has the same sign as its electrostatic counterpart 
(negative and attractive), while the effective triplet potential has the oppo- 
site sign (positive and repulsive). The induced pair attraction originates from 
the depletion of b particles from the space between pairs of closely approach- 
ing a particles and the resulting imbalance in osmotic pressure. Furthermore, 
in contrast to systems of charged particles with electrostatic interactions, in 
systems of hard particles with excluded-volume interactions, the volume term 
depends only trivially on the density of the a particles and the effective many- 
body interactions do not generate corrections to lower-order effective inter- 
actions. One must take care, however, not to draw general conclusions, since 



(116) 



(117) 
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the neglected higher-order terms in (113) can significantly modify the effec- 
tive interactions - especially in concentrated systems - introducing density- 
dependence and even changing the sign. In binary hard-sphere mixtures, for 
example, packing of smaller spheres around larger spheres induces effective 
pair interactions between the larger spheres that can exhibit a repulsive bar- 
rier and even long-range oscillations [39, 69]. 

A simple, but important, example of a binary mixture of spheres is the 
Asakura-Oosawa-Vrij (AOV) model [70, 71] of mixtures of colloids and free 
(nonadsorbing) polymers. The AOV model treats the colloids (a ^ c) as hard 
spheres, interacting via an additive hard-sphere pair potential Vcc{r) (HI), 
and the polymers {b ^ p) as effective, coarse-grained spheres that have hard 
interactions with the colloids, 

''''P^'^' \0, r>Rc + Rp, ^ ' 

but are mutually noninteracting (ideal): Vpp{r) = for all r. The radius Rp 
of the effective polymer spheres is most naturally identified with the polymer 
radius of gyration. The neglect of polymer-polymer interactions is strictly 
valid only for theta solvents [11], wherein monomer- monomer excluded- volume 
interactions effectively vanish. In this special case, the polymer in the system 
behaves as an ideal gas confined to the free volume, Vf = aV, where the 
free volume fraction a is defined as the ratio of the volume available to the 
polymer centres (i.e., not excluded by the hard colloids) to the total volume. 
At equilibrium, equality of the polymer chemical potentials in the system, 
IJ-p = ksTln^npAp/a), and in the reservoir, ii^^ = kBTlii{n^p^ Ap), implies 

(r) 

that the corresponding polymer densities must be related via Up = arip . 
This simple relation imposes equality also of the polymer osmotic pressures 
in the system and reservoir, Up = Hp ' = ksTrip ', thus rendering (114) and 
(117) exact. 

Within the AOV model, not only are the effective interactions exact, but 
the effective one-body and pair interactions have analytical forms, since the 
single-sphere excluded volume is simply V^xc = (47r/3)(i?a -|- Rb)^, and the 
convex-lens-shaped pair intersection has a volume 

FX')(r) = |f (^^' + ^^) 

[ 0, r>ac + dp, 

(120) 

where = 2i?c and ap = 2Rp are the particle diameters and q = Rp/Rc 
is the size ratio. The effective one-body interaction (116), being linear in the 
colloid density, does not affect phase behaviour, but does contribute to the 
total osmotic pressure. The effective pair potential described by (117) and 
(120) consists of a repulsive hard-sphere core and an attractive well with 
a range equal to the sum of the particle diameters and depth proportional 



3r 



2((Tc + (Tp) ^ 2(c7e + ap) 



ac <r < ac + (7p 
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to the reservoir polymer osmotic pressure. For sufficiently small size ratios 
(q < 0.154), such that the spherical exclusion spheres surrounding three col- 
loids never intersect, the effective triplet (and higher-order) interactions are 
identically zero. For q > 0.154, the effective triplet interaction is nonzero and 
repulsive. 

Although limited to ideal polymers, the AOV model gives qualitative in- 
sight even for interacting polymers. For sufficiently large size ratios {q > 0.45) 
and high polymer concentrations, polymcr-dcplction-induccd cfFectivc pair at- 
tractions can drive demixing into colloid-rich (polymer-poor) and colloid-poor 
(polymer-rich) fluid phases [8]. Simulations of the effective Hamiltonian sys- 
tem [38] and of the full binary AO model [72] indicate that for smaller size 
ratios (g < 0.45), fluid-fluid demixing is only metastable, being preempted 
by the fluid-solid (freezing) transition. With the phase diagram of the AOV 
model now well understood, recent attention has turned to exploring the phase 
behaviour of mixtures of colloids and interacting polymers (Sec. 4.3). 

Cluster Expansion Approach 

An alternative, and elegant, approach to modelling effective interactions in 

systems of hard particles has been developed recently by Dijkstra et al [38, 39]. 
This powerful statistical mechanical method is similarly based on integrating 
out the degrees of freedom of one species of particle in the presence of flxed 
particles of another species. The essence of the method is perhaps most trans- 
parent in the context of the AOV model of colloid-polymer mixtures [38], 
deflned by the Hamiltonian H = K + Hcc + Hep, where K is the kinetic 
energy and 



describe the colloid-colloid and colloid-polymer interactions. Assuming an in- 
finite reservoir that exchanges polymer with the system, the natural choice of 
ensemble is the semigrand canonical ensemble, in which the number of colloids 
Nc, volume V, temperature T, and polymer chemical potential fj,p are fixed. 
Within this ensemble, which treats the colloids canonically and the polymers 
grand canonically, the appropriate thermodynamic potential is the semigrand 
potential f2, given by: 



polymer coordinates, Zp = e^^p/A^ is the polymer fugacity, and and A, 




(121) 



exp(-/3r2) = {{exp[-P{Hcc + Hcp)])p)c 
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are the colloid and polymer thermal wavelengths. Here i?eff = H^c + is the 

effective one-component Hamiltonian, where Q.p represents the grand potential 
of the polymers in the presence of fixed colloids, defined by 



exp(-/3f2p) = (exp(-/3i?cp))j 



EZp_ 

JV„=0 



-Op 



dr exp 



-P^Vcp{\T^ 
i=l 



exp i-^p J 



i=l 



(123) 



Equating arguments of the exponential functions on the left and right sides, 
we have 



-I3n, 



dr exp 



■r|) 



(124) 



As shown by Dijkstra et al [38, 39], the polymer grand potential can be 
systematically approximated by a cluster expansion technique drawn from the 
theory of simple liquids [20]. Defining the Mayer functions 



fij = exp[-/3ucp(ry)] " ^ = | 



-1, 
0, 



Tij <C Rq Rp 
^ij ^ R-c Rpi 



(125) 



(124) can be expanded as follows: 



= zp [dr J + 

= zp h + E ^^i^^o + ■■■] 

^ \ i=l i<fe=l / 



(126) 



Equation (126) is a form of cluster expansion well-suited to systematic approx- 
imation by diagrammatic techniques [20] . Successive terms in the summation, 
generated by increasing numbers of colloids interacting with the polymers, 
correspond directly to the one-body volume term and effective pair and many- 
body interactions. In fact, from (126), the semigrand potential can be written 
in a form that is precisely analogous to the free energy defined in (114): 

/ . Na Na 

f3Qp = m+zp AT.yW V^^Hr.j) + ^ E V^^Hrij,rik) - ■ 

(127) 

where /3f2o = —ZpV is the grand potential of the reference system of pure 
polymer. In the case of ideal polymer, the polymer fugacity is simply related 
to the reservoir osmotic pressure via Zp = (3n^\ Note that (114) and (127) 
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are completely equivalent, differing only with respect to the relevant ensemble, 
with (114) applying in the canonical ensemble and (127) in the semigrand 
canonical ensemble. 

For mixtures of colloids and polydisperse polymers, (126) generalizes to 

= ^24"^ I dr, n(l + f\^), (128) 

fe ^ i=l 

where Zp'^'' is the fugacity of polymer species k and f^^^ is the corresponding 
Mayer function. From (128), the effective pair potential is then simply a sum 
of depletion potentials, each induced by a polymer species of a different size. 
On the other hand, in the case of interacting (nonidcal) polymers [vpp{r) ^ 0], 
the effective interactions - even for monodisperse polymers - are considerably 
more complex. Although (123) then can be formally generalized to 

exp(-/3J7p) = (exp[-/3(iJ,p + ifpp)])p 

= E ^ / dr"^ n 11(1 + 11(1 + /If')' (129) 

JVp=0 P' •' i=lj=l k<l 

where f^'^^ and /fj^ are the Mayer functions for colloids and polymers, respec- 
tively, practical expressions are less forthcoming. Using diagrammatic tech- 
niques, Dijkstra et al [39] have further analysed (129) and demonstrated its 
application to the phase behaviour of binary hard-sphere mixtures. 

4.3 Systems of "Soft" Particles 

Soft particles arc here defined as macromolcciilcs having internal conforma- 
tional degrees of freedom. Prime examples are flexible polymer chains (linear 
or branched), whose multiple joints allow for many distinct conformations. 
While bare monomer-monomer interactions - cither intrachain or interchain 
- can be modelled by simple combinations of excluded- volume, van der Waals, 
and Coulomb pair interactions, the total interactions between long, fluctuat- 
ing chains can be highly complex, rendering explicit molecular simulations of 
polymer solutions computationally challenging. 

Among several practical techniques developed for averaging over the inter- 
nal structure of soft particles to derive effective pair interactions, we discuss 
here two that are speciflcally suited to polymers in good solvents. A more 
extensive survey is given in the review by Likos [28] . One method is based on 
the general principles of polymer scaling theory [11], which describes the prop- 
erties of polymers in the limit of infinite chain length, i.e., segment number 
N ^ oo. The starting point is a formal expression for the effective pair inter- 
action between the centres of mass, at positions Ri and R2, of two isolated 
chains (labelled 1 and 2): 
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/3veff(i?i2) = -ln(^^ J Vvi j Pr2Pcm(Ri)Pcm(R2)exp[-/3ff[{ri},{r2}]^ , 

(130) 

where Zi is the partition function of a single isolated chain, J Dr^ represents 

a functional integral over all conformational degrees of freedom of chain a, 
Pcm(R-a) denotes the number density of the centre of mass of chain a, and 
the Hamiltonian iJ is a functional of the conformations, {ri} and {r2}, of the 
two chains. In the special case of two chains, with one end of each chain fixed 
and the two fixed ends separated by a distance r, (130) becomes 

where Z2{r) is the partition function of the constrained two-polymer system 
and £^2(00) is the same in the limit of infinite separation. Scaling arguments 
[73] suggest that in the limit N 00. in which the only relevant length 
scales are the separation distance r and the polymer radius of gyration 
Z2(r)/Z2(cx)) (X {r/RgY , where a; is a universal exponent. It follows that 

I3v,s{r) oc - In (^^^ , r < Kg, (132) 

which describes a gently repulsive effective pair potential. 

Similar scaling arguments have; been extended to star-branched polymers 
[28, 73], consisting of linear polymer chains (arms) all joined at one end to 
a common core. The same form of effective pair potential results, but with 
an amplitude that depends on the number of arms and the solvent quality. 
A refined analysis by Likos et al [74] leads to an explicit expression for the 
effective pair potential between star polymers, which extends (132) to r > Rg, 
specifics the prefactor, and is consistent with experimentally measured static 
structure factors. 

An alternative coarse-graining approach, well-suited to dilute and semidi- 
lute solutions of polymer coils in good solvents, is based on the physically 
intuitive view of polymer coils as "soft colloids" [75], whose effective interac- 
tions may be approximated using integral-equation methods from the theory 
of simple liquids [20] . The basis of this approach is the fmidamental relation 
between the pair distribution function g{r), direct correlation function c(r), 
and pair potential v(r) of a simple liquid: 

g{r) = exp[-Pv{r) + g{r) - c{r) - 1 - &(r)]. (133) 

Equation (133) follows directly from the Euler-Lagrange relation (63) for the 
nonuniform density of particles ng{r) around a central particle, where the pair 
potential plays the role of the external potential and the bridge function b{r) 
subsumes all multiparticle correlation terms. Practical implementation begins 
with a numerical calculation of g{r) between the polymer centres of mass, 
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e.g., by molecular simulation of explicit self-avoiding random- walk chains, and 
proceeds through inversion of (133) to determine the effective centre-to-centre 
pair potential v{r). 

The inversion of g{r) combines (133) with the Ornstein-Zernike relation 
(64) and an approximate closure relation for the bridge function. Louis and 
Bolhuis et al [75, 76] have demonstrated that the HNC closure, b{r) = 0, 
gives an accurate approximation for the effective interactions. The resulting 
softly repulsive, Gaussian-like, effective pair potential has a range comparable 
to the polymer radius of gyration, an amplitude ~ 2A:bT, and is only weakly 
dependent on concentration, implying the relative weakness of effective many- 
body interactions [77]. In self-consistency checks, simulations of simple liquids 
interacting via the potential v{r) are found to reproduce, to within statistical 
errors, the same centre-to-centre g{r) as simulations of explicit chains. The 
same method has been applied also to calculate the effective depletion-induced 
interaction between hard walls [75, 76] and between colloids in colloid-polymer 
mixtures [78, 79]. 

The scaling and soft-colloids approaches both reach the common conclu- 
sion that effective pair interactions between polymers in good solvents are 
ultrasoftly repulsive. As the centres of mass of two polymers approach com- 
plete overlap, the pair potential between star polymers diverges very slowly 
(logarithmically), while that between linear chains actually remains finite. 
These characteristically soft effective interactions contrast sharply with the 
steeply repulsive short-ranged pair interactions between colloidal particles, 
going far in explaining the unique structural properties and phase behaviour 
of polymer solutions. 

5 Summary and Outlook 

The key message of this chapter is that soft materials, comprising complex 
mixtures of mesoscopic macromolecules and other microscopic constituents, 
often can be efficiently modelled by preaveraging over some of the degrees of 
freedom to map the multicomponent mixture onto an effective model, with 
fewer components, governed by effective interparticle interactions. In gen- 
eral, the effective interactions are many-body in nature and dependent on the 
thermodynamic state of the system. Briefly surveyed were several recently 
developed statistical mechanical methods, including response theory, density- 
functional theory, and distribution function theory. These powerful methods 
provide systematic and mutually consistent approaches to approximating ef- 
fective interactions, and have broad relevance to a variety of materials. Specific 
applications were illustrated for electrostatic interactions in charged colloids 
and excluded-volume interactions in colloid-polymer mixtures. 

Effective interactions are often simply necessitated by the computational 
impasse presented by fully explicit models of complex systems, especially soft 
matter systems with large size and charge asymmetries. At the same time, 
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however, effective models can provide conceptual insight that may be diffi- 
cult or impossible to extract from explicit models. Consider, for example, the 
subtle interplay of entropy and electrostatic energy in charged colloids or the 
important role of polymer depletion in colloid-polymer mixtures, effects that 
arc elegantly and efficiently captured in effective interaction models. While 
computational capacity will likely continue to grow exponentially in coming 
years, conceptual understanding of soft materials will also continue to benefit 
from the theoretical framework of effective interactions. 
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